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I. INTRODUCTION 



Fuel tanks comprise the greatest portion of vulnerable target 
volume in modern military aircraft. This degree of vulnerability of the 
fuel tanks is of paramount importance in determining aircraft surviva- 
bility. Fuel tanks which are subjected to ballistic impact can sustain 
severe damage as a result of the phenomenon known as hydraulic ram. For 
our purposes, hydraulic ram is defined as the production of pressure 
loadings by ballistic impact and their effects on the fuel tank and 
its components. The term hydraulic ram is a misnomer in that it sug- 
gests a quasi- static compression process. Actually the mechanism of 
hydraulic ram consists of several phases. The shock phase is generated 
when the projectile penetrates the tank wall. Projectile penetration 
produces a stress riser in the tank wall from which cracks can propa- 
gate radially. Shock wave formation in the fuel due to projectile 
penetration into the fluid causes very high local pressures which may 
be sufficient to cause failure of the entry wall in the neighborhood of 
the entry point. As the projectile moves through the fluid it leaves a 
vapor- filled void caused by flow separation from the projectile surface. 
During this so-called cavity phase the projectile's kinetic energy is 
transferred to the fluid by cavity formation. Cavity growth is subse- 
quently defeated by hydrostatic pressure and the cavity collapses. 

The collapse is not total because of fluid vapor compression. The 
cavity oscillates until equilibrium is reached in the fluid. Cavity 
oscillations produce corresponding pressure loadings on the tank walls 
which can cause catastropic failure of tank components. The exit wall 
of the tank is also heavily loaded as the projectile passes through. 

The purpose of the investigation to date has been to develop an 
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understanding of the hydraulic ram phenomenon. Fuel tank response to 
ballistic penetration cannot be fully understood without comprehension 
of the individual hydraulic ram components and their interaction. This 
understanding can be obtained by conducting experiments and developing 
analytical- numerical tools which predict the events that occur during 
the damaging processes. With this understanding and the analytical- 
numerical tools available, fuel tanks which minimize hydraulic ram can 
be designed. Once fuel tanks are protected against this type of threat, 
aircraft vulnerability will be lessened. 
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II. ANALYTICAL AND NUMERICAL STUDIES 



A. A Study of the Wall Motion, Treated as a Single Degree of Freedom 
System (Development of Piston Theory) 

Consider a typical portion of a fuel tank wall treated as a single 
degree of freedom system as shown in Fig. II-A-1. The area of the por- 
tion of the wall is A and the total mass is M. The resistance to 
displacement provided by the surrounding wall is modeled by the spring, 
whose constant is K. The equation governing the elastic motion of the 
wall is assumed to be 



Mw + Kw = Ap( t) (II-A-1) 

where p(t) is the total time dependent pressure on the wall due to the 
incident pressure waves in the fluid and to the waves reflected and 
radiated from the wall and w is the wall displacement. The dots indi- 
cate derivatives with respect to time t. Assuming that the fluid 
satisfies the linear wave equation and applying a step pressure pulse 
normal to the wall as the incident wave gives 



p(t) =0 t < 0 
p(t) = P o - t s 0 



( II-A-2 ) 



where p Q is the pressure in the incident wave, p is the fluid density, 
and P is the velocity potential of the reflected and radiated waves. 
The subscript i denotes the value of p at the wall. The potential P 
is governed by the wave equation 

c P = P (II-A-3) 
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where c is the acoustic velocity in the fluid and primes denote deriva- 
tives with respect to x, the axial space coordinate normal to the wall. 

\ 

At the wall the velocity of the fluid is equal to the velocity of the 
wall . Thus , 

w = v. + ft ' (ii-a-4) 



where v. is the velocity of the incident wave and is equal to p Q /(pc). 

Equations (l), (2), (3), and (4) govern the behavior of the fluid- 
structure interaction process. The solution to Equation (3) for p can 
be given in the form 



P = ae i^.(x+ct)/a + ge iX(x-ct)/a 



(H-A-5) 



where a and 0 are arbitrary constants, X is a scalar, and a is the 
value of x at the wall. The first quantity on the r.h.s. of Equation (5) 
represents waves traveling away from the wall, and the second quantity 
represents waves traveling toward the wall. Note that for the waves 
traveling away from the wall, i.e., the reflected and radiated waves, 

# I 

p = cp (ii-a-6) 



Thus, 




according to Equation (4). Hence 

f> = c(w-v ) (II-A-7) 

Substituting p given by Equation (7) into Equation (2) gives 
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(ii-a-8) 



p(t) = P o + pc^-w) t > 0 



Examination of Equations (l) and (8) reveals that the variables fi and 
w have now been separated, i.e., the response of the wall can be compu- 
ted based upon the pressure and velocity of the incident pressure wave. 
Application of Equation (8) to the more general two-dimensional wall 
response is known as the piston theory (Refs. 1 and 2). Substituting 
Equation (8) into Equation (l) leads to 




(II-A-9) 



where 







The solution to Equation (9) can be given in the form 




(lI-A-10a) 



w = w 



st 




e 



-C<rt 




(11-A-lOb) 
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(II- A- 10c) 
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w 



Ap p 
*0 _ o 

st " K 2 

mao 



when the wall is at rest prior to the application of the step pressure. 

If the contribution of the spring is neglected, we have a one- 
dimensional pressure wave reflecting off of an infinite wall. Hence, 

1 pc 

ud = 0, £ = », ^ = 2 m Problem was considered in Ref. 3, where 

the equation for the pressure on the wall was taken as 



P(t) = P o e 



pet 

m 



(Il-A-I0d) 

Typical values for c, p and m for hydraulic ram in fuel tanks are 

c = 53,600 in./sec 

p = 7*64 x 10 ^ #-sec^/in.^ 

5 2 ^ 

m = 2.59 x 10 #sec / in. (0.1 in. aluminum wall) 

= 3.64 x 10 ^ #sec^/in.^ (0.5 in. steel wall) 



Therefore 



— = 0.158 x 10^ (0.1 in. aluminum wall) 

m sec 



— = 0.01125 x 10 6 — (0.5 in. steel wall) 

m sec 

Consequently, when uj < 15000 rad/sec for the 0.1 in. aluminum wall 
and uo < 1000 rad/ sec for the 0.5 in. steel wall, £ > > 1, i.e., the 
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system is heavily damped, and Equations (10a - 10c) simplify to 



w = 






UUt 

~Q 



uut 



2C 




w 



w 



st 




uut \ 

e 2C e -2C^t 



(II- A- 11a) 



(II- A- lib) 



P = P„ \ 1 - e 



uut uut 
2 C + 2Q e -2C«*| 



( II-A-llc) 



A plot of — versus time is given in Fig. II-A-2a for — = 0.2, 0.1, 
p ° 6 l 

0.05, and 0.01 x 10 , and 1 usee. In this figure the assumption is 

o sec 

d. d ~\ 

made that < 100, i.e., uu < 1,000 rad/sec for — = 0.01 x 10 

pc 5 5 ’ ' m sec 

and uu < 4,^70 for — = 0.2 x 10^ — — . Under these conditions, Equation 
m sec 

(lie) for p simplifies to Equation (lOd), provided t < 500 psec. Curves 
of this type are also given in Ref. 3> Fig. 8, for a 1 inch plexiglas 
wall and a 0.063 inch aluminum wall with an exponentially decaying inci- 
dent pressure wave. 

For t > 500 psec the elastic property of the wall begins to become 

effective and reduces the wall velocity. This causes the pressure to 

2 

build up. A plot of is given in Fig. II- A- 2b for = 100, 20, 10 

± P o pC 

and 5 ? and 100 psec < t ^ 100,000 psec. For these conditions, 

sec 



Equations ( 11a - 11c) simplify to 

'st V 



w = w at 1-e 






uut / ujt 

2C _ f. 2C, — 2 pout ] 

, P = P Q ll-e e K 1 



(II- A- 12) 



Examination of Figs. II-A-2a and 2b reveals that for the values of 
the parameters considered, the inertial effects and the elastic effects 
occur in essentially separate periods of time. This is due to the 
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assumption on uo. When the stress wave strikes the wall, the wall de- 



flects as if there were no elastic effects. After approximately 500- 
1000 ^sec the wall starts to slow down due to the elastic effects, the 
pressure builds up, and eventually the wall comes to rest. The motion 
of the wall is that of a heavily damped oscillator* . This phenomenon 
holds provided 



2 

moo 

m 



< 100 



1 

sec 



For 



— = 0.01 x 10^ — — , uo < 1000 rad/sec 

m sec 5 ’ 



The lowest natural frequency of a square plate of width L is 



(JO = 



2tt 



V 12(l-» 2 )p 



(II- A- 13) 



where E is Young's modulus, h is the wall thickness, and u is Poisson's 
ratio. For both aluminum and steel walls 



0) 



1.2 x 10 ^ (—x\ rad/ sec 



(lI-A-14) 



For 



L = 20 in. and h = 0.5 in. 



U) 



150 rad/ sec 



For 



L = 10 in. and h = 0.5 in. , u> 600 rad/ sec 



*The natural frequency for critical damping is 79*000 rad/sec (12,600 
Hz) for the 0.1 aluminum wall, and is 5,620 rad/sec (884 Hz) for the 
0.5 in. steel wall. 
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Thus, uo < 1000 rad/ sec is a realistic limit for the fundamental frequency 
of fuel tank walls. 

Since the elastic response takes place relatively late in the inter- 
action process, the question arises as to the possible validity of 
treating the stress wave as an impulse loading. For example, suppose 
the incident pressure pulse at the wall is of the type shown below 




The solution for the wall displacement (provided £ > l) for t > t is 



w 




e 



( t— T )uj 

2 £ 





e 



cu( t- t) 

2 £ 



- 2 £u>( t- t ) 
(II- A- 15) 



t > T 



where w^ and w^ are the displacement and velocity of the wall at t = t. 
The first term in Equation (15) is the long time response; the second 
term damps out very quickly. If the stress wave pulse is treated as an 
impulsive loading, the displacement and velocity at t = t are 



w = 0 

T 




Hence 



w impulse 



V 

2£um 




U>( t- T ) 
2 £ 



• U)( t- T ) 

2 £ 

e 



■2£u>( t-T) 



(II-A-16) 
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The correct solution for w using Equations (10a) and (10b) to define 
w and w is 

T T 



w 



T 





w 






and 




mu) 



(JUT . 



1-e 



U}( t-T) 




2Q 

e * e 



HI . zim \ i%l) -2;o,(t-T)| 

(II- A- 17) 



t > T 



provided Q > > 1. When t is such that 



UUT 

2C 



ojt 

' 2G 



i.e., t < < , Equation (17) simplifies to 



muo 



w 



2£um 



uu( t- t ) 

2G 



-L -2Gwt -2£uu(t-T)^ 



+ — e 

H 2 



2Q 

e e 



t > T 

(II-A-18) 



Thus, comparing Equations (l6) and (l8) reveals that the long time 
response is correctly approximated by the impulse assumption, but the 
short time response is inaccurate. 

The conclusions that can be drawn from this analysis are: 

1. The wall responds to the pressure wave as a very heavily 

damped oscillator provided the natural frequency of the wall 
is less than 1000 rad/sec (rEL50 Hz). 
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2. The pressure on the wall drops off to a small percentage of 
the incident pressure within the first 100 usee. It begins to 
build up when t > 1000 usee. Thus, the reflected pressure from 
the wall is essentially the same as the reflected pressure from 
a free surface for the period 100 usee < t < 1000 usee. 

3. If the incident pressure is on for a finite length of time t 

such that t < < — ( — = 0.158 x 10^ — — for a 0.1 in. aluminum 
pc \m sec 

wall, — = 0.01125 x 10^ — — for a 0.5 in. steel wall) the long 
m sec / 

time response is accurately predicted but the short time response 
is inaccurate. 
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Figure II-A-2b Wall Pressure Vs. Time (x > 100 y sec. 



B. Entry Wall Response to the Shock Phase 



A preliminary study of the structural response of the entrance wall 
of a fuel tank during the initial shock phase of the hydraulic ram has 
been conducted. The solution for the pressure in the fluid derived by 
Yurkovitch (Ref. 4 ) was reduced to a simple analytical form for use in 
the structural analysis digital computer program SATANS*. The form 
selected for the pressure p at a radius r and time t was 




where 1 > { } > 0. This form corresponds to a peak pressure that varies 
t 

as p* — with a linear decay behind the front to zero pressure. The 
t l 



location of the front is given by C r 



ft)’ 



For a 0.50 caliber project- 



ile impacting H^O (Fig. 22 of Yurkovich's report) 



p* = 150 ksi 

t^ = tg = t^ = 10 psec 

C = 1.22 

C 2 = 1.7 in. 

m = 1.2 



*SATANS computes the transient, geometrically nonlinear response of 
elastic, axisymmetric plates and shells subjected to arbitrary loads. 
It is described in detail in Refs. 5 and 6 . 
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n = 0.87 
q = 0.68 

The pressure loading given by this model is shown in Figure II-B-1 and 
was applied to a circular aluminum plate with a radius of 10 in. and 
a thickness of 0.1 in. The outer edge of the plate was assumed to be 
clamped. The results for the plate displacements computed by SATANS 
are shown in Fig. II- B- 2 for t = 5? 10, 20, 30, and 50 usee. The 
maximum membrane force in the plate was approximately 60,000 lb/in. 
and occurred at the center of the plate at t sa 35 Usee. The maximum 
stress caused by the force and by the bending moment was 600,000 psi 
and 720,000 psi respectively. 

The analysis for the response of the entrance wall during the shock 
phase discussed above neglects several important features. First, it 
does not take into consideration the fact that the motion of the wall 
relieves, to some extent, the pressure acting on the wall. Second, 
the very large stresses developed in the wall are considerably higher 
than the yield stress and the rupture stress of aluminum. Hence, 
yielding and cracking will occur. Third, the effects of shear defor- 
mations and rotary inertia have not been considered. Of these three 
features, the effect of wall motion upon the loading pressure is consi- 
dered to be of prime importance in the initial phase of this program. 

An exact method for determining the total fluid pressure within a 
closed, flexible tank has been developed by Lundstrom (Ref. 6), This 
model was considered as the algorithm for the pressure loading. 

However, due to difficulties encountered in the programming of the 
model, Lunstrom recommended that the piston theory be used in place 
of his exact model. Piston theory is described in II-A. Application 
of this theory uncouples the structure and fluid equations and allows 
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the damped response of the wall to be computed using the incident 
pressure as the load on the wall. This approximation has been in use 
since the early 1940' s when it was applied to the study of the effect 
of underwater explosions on ship's plates. 

The theory is based upon the assumption that the normal pressure on 
the wall is given by Equation II-A-8, i.e., 

p(t) = p + p (v. - w) 
r o c i 

which is the equation for one- dimensional fluid flow. Since SATANS did 
not contain the damping term pew, it was necessary to modify the code 
to incorporate this procedure. 

The results for the response of the same circular plate subjected 
to the same pressure pulse described above but now with piston theory 
added, are presented in Fig. II-B-3 at several times. Comparing the 
two figures, II- B- 2 and 3, reveals that accounting for fluid- wall inter- 
action has a considerable effect on the wall response. Nevertheless, 
for this example the wall stresses are still above the yield, i.e., 
the maximum membrane force is 8,900 lb/in. and occurs at the center at 
t = 15 psec, the maximum bending moment is 478 in. -lb/in. and occurs at 
r = 1 in. at t = 15 i^sec, the maximum stress is 32? ksi and occurs at 
r = 1 in. at t = 15 qsec. 
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Pressure - p x 10 (psi) 





FIGURE E-B-l 



SHOCK PRESSURE DISTRIBUTION FOR 
.50 CALIBER PROJECTILE IMPACTING 
H 2 0 
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Displacement ( in.) 




FIGURE H-B-2 NORMAL DISPLACEMENT vs. RADIAL 

DISTANCE WITHOUT FLUID INTERACTION 
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Displacement (in.) 




FIGURE IT-B-3 DISPLACEMENT OF FUEL TANK 

WALL , WITH PISTON THEORY 
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C. Modification of the Northrop Internal Blast Code (BR-l) to Account 



for Hydraulic Ram Loading and Fluid- Structure Interaction 

The Northrop Corporation is developing a finite element digital 
computer code for predicting the transient, inelastic, large deflection 
response of flight vehicle structures subjected to internal blast loading 
and fragment penetration. Ball has been following the development of 
this code with the objective of modifying it to account for the hydraulic 
ram loading and fluid- structure interaction that occurs when bullets 
and metal fragments penetrate aircraft fuel tank walls. The main features 
of the program are summarized below: 

1. The elements are the isotropic flat rectangular plate and the 
straight beam. The cross-section of the beam element is either rectan- 
gular, I, or T-shaped. 

2. The Bernoulli- Euler hypothesis is used for both the plate and 
the beam element. 

3. The nodes of the plate element are at the four corners. The 
degrees of freedom at each node are the three displacements u, v, and w 
and two rotations 6^ and 9^. A linear edge displacement state is assumed 
for the inplane displacements u and v, and a cubic function is assumed 
for the normal displacement w. 

4. Large deflection strain- displacement relationships are used for 
the plate element that are valid for strains as large as 10 % and rota- 
tions as large as 20°. 

5. Gaussian quadrature is used for the integration process. Three 
points in each direction are used. 

6. The displacement state assumed for the beam element leads to an 
axial force, two transverse shearing forces, and a uniform torsion over 
the length of the element. 
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7. Warping of the beam due to torsion is neglected. 

8. The constituitive relations for the plate and beam material 
incorporate the von Mises yield condition, the Prandtl-Reuss flow rule 
and isotropic work hardening. The material behavior is specified by a 
piecewise linear uniaxial stress-strain relationship. The plate cross- 
section is divided into four layers for the stress computations. 

9. Offset points are used to relate beam element nodal to plate 
element nodal deflections. 

10. An initial state of stress feature is available wherein each 
element will be in a uniform state of stress with all deflections zero. 

11. Permissible loadings are distributed (pressure) loading normal 
to the deformed structure, either continuous or impulsive, impulsive 
concentrated forces at the node points, and an automated pressure frag- 
mentation impulse loading from the NOL BLAST Code. 

12. Account is taken of randomly located small perforations of a 
plate element due to fragment penetration by reducing the membrane and 
bending stiffnesses, the pressure load, and the mass. 

13. The timewise numerical integration scheme is an explicit proce- 
dure . 

14. The lumped mass technique is used to develop the mass matrix. 
Rotary inertia is considered. 

15. An eigenvalue routine is incorporated to automatically deter- 
mine the numerically stable time increment. 

16. The program has a restart capability. 

17. The code is operational on the CDC 6600 (FORTRAN EXTENDED) sys- 
tem. A structure with a maximum of 300 joints and 400 members can be 
analyzed. A total core storage of 225,000 words OCTAL (CDC system) or 
270,000 bytes (IBM 320-165 system) is required. 
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The modifications to the BR-1 Code to make it applicable to bullet 
and fragment penetration into fuel tanks can be divided into two main 
categories: (l) incorporation of the hydraulic ram pressure loading, 

and (2) accounting for the fluid- structure interaction phenomena. Due 
to the difficulties involved in treating the fully coupled problem, the 
fluid- structure interaction phenomena will be approximated by the piston 
theory described in II-A. Thus, since the BR-1 Code does not account 
for damping, it will be necessary to modify the code by adding the 
damping term pew to the equations of motion. 

In the piston theory, the pressure on the structure is p Q + pev^ 
where p^ is the pressure in the incident wave, p is the fluid density, 
c is the acoustic velocity in the fluid, and v^ is the fluid velocity. 
Lundstrom has developed a computer code that will predict the pressure 
and velocity throughout a rectangular volume of an irrotational, ideal 
compressible fluid due to a penetrating bullet. The boundaries of the 
fluid may be either free or rigid. This code will be used to determine 
the pressure on the fuel tank as input data to the BR-1 code. 
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D. Analysis of Cavity Motion 



In order to gain insight to forces and velocities associated with 
cavity motion two problems were formulated and solved by Holm (Ref. ,7). 
One was the collapse of a spherical cavity, an event which occurs when 
a vacuum exists in the cavity. The other was the oscillation of a 
spherical cavity with entrapped gases. 



Cavity Collapse 

The equations of motion for an incompressible fluid or density p can 
be manipulated into the form 



R 




2 

R d R 



(iI-D-l) 



where R is the cavity radius. The pressure within the cavity is p, and 
at infinity p^. At time t = 0 the radius is R q . For the cavity collapse 
problem p = 0, and equation (l) can be manipulated into the following 
form 




(II-D-2) 



where X = R/R 0 * Note that the factor multiplying the integral has the 
dimensions of time and defines a characteristic time 



T 




(II-D-3) 



Two examples are shown in Figures 1 and 2. For both cases the collapse 
time is t = 0.6 t. in terms of the characteristic time, the curve in 
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either figure is a universal curve for this problem. As a rule-of-thumb, 
cavity collapse is about 60 % of the characteristic time. 



Cavity Oscillation 

When p / 0 in Equation (l), there occurs oscillations of the cavity. 
It is assumed that the compression and expansion of the trapped gas is 
adiabatic and lossless, i.e., isentropic. Consequently, pressure and 
cavity radius are related by 



P 



P 



o 




3\ 



(ii-d-4) 



where p Q is the gas pressure when R = R q . The ratio of heat capacities 
0^ /C v is symbolized by y. Introduction of Equation (4) into Equation 
(l) yields 
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T 




(II-D-5) 



where a new parameter a has been introduced. This parameter depends on 
initial ratio of pressures within the cavity and at infinity and is 
defined by 



a 




(II-D-6) 



In Equation (5) Z equals R q /r which should be noted is the reciprocal 
of X in Equation (2). A solution of Equation (5) is shown in Fig. 3. 
The complete curve cannot be obtained from a single integration. When 

3 3y 

Z (cv + l) - a Z - 1 = 0, the integrand approaches infinity. This 
occurs when a maximum or minimum in radius or pressure is attained. 
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At this point a new set of R and p values are introduced and inte- 

o o 

gration proceeds. 

For the example of Fig. 3 , the pressure starts at p^ /p^ = 6.67 and 
Z = 1. When p/p^ = 0.4 at Z = .80 the integrand approaches infinity. 
This is at the minimum of the p/p^ trace. One continues at this point 
letting p /p =0.4 and a new Z defined as 1.000. Integration then 

o ' oo 

continues to get the portion of the curve from minimum p/p^ to maximum 
p/p^, i.e., the part of the curve with positive slope. For the case 
of a - 16.675 5 one cycle of oscillation requires a time of 1.35 T . 

In the interpretation of Fig. 3 note that large p/p^ occurs for 
small R or large Z. Although 0.4 £ p/p ^ 6.67 in Fig. 3* the range of 
radius is 0.64 £ R/R q - 1.24. Large variations in p/p^ are caused by 
small changes in R/R^ due to 3y power dependence. 
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Time milliseconds 

Fig. I-D-l Cavity Collapse Time from an 
Initial Radius of 0.5 Foot. 
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Fig. IL-D-2 Cavity Collapse Time from an 
Initial Radius of 1.0 Foot. 
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t/T 

Fig. n-D-3 Hydraulic Ram Cavity Pressure 
Oscillation. 
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III. HYDRAULIC RAM EXPERIMENTAL FACILITIES 



A. Introduction 

The design, fabrication, and calibration of a ballistic range for 
aircraft survivability studies was accomplished under the present con- 
tract. The range was designed specifically with the hydraulic ram 
problem in mind. Of paramount importance was a consistently accurate 
method for determining a time history of the projectile velocity before 
projectile-tank impact and during hydraulic ram events that occur after 
penetration. Since projectile attitude at entry is important, shadow- 
graph stations were required along the flight path to make an estimate 
of projectile attitude. Two shadowgraph stations were also required to 
picture the hydraulic ram phenomenum as a function of time within the 
tank. The shock wave system, projectile position, and cavity shape and 
size could be determined using this data. The electronics associated 
with the ballistic range was designed to give stable, repeatable measure- 
ments at very small time intervals. Calibration of the present range 
components has shown that consistently accurate data is obtained. As 
a result, a reliable means of obtaining projectile energy as a function 
of time during hydraulic ram has been accomplished. 
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B. Ballistic Range Components 



Basic elements of the ballistic range are shown in Figure III-B-1. 
The rifle mounting system is composed of the rifle mount and stand. The 
rifle mount was adjustable in azimuth ( + 8 degrees) and elevation 
(+3 degrees to - 5 degrees) for ease in boresighting. Rifles of 22.2 
and 30 caliber may be used. Figure III-B-2 shows the complete rifle 
mounting system and Figure III-B-3 shows the rifle mount in detail. 

The shadowgraph stations are composed of a bullet sensor , time delay 
unit, spark source, collimating lens, and a shadow box with reference 
grid for mounting the Polaroid film holder. Figure III-B-4 shows a 
typical shadowgraph station. The bullet sensor is a chronograph screen 
that has a five volt D.C. signal shorted to ground across it. The 
bullet sensors provide start and stop pulses to the counters and delay 
circuits. The delay units generate time delayed pulses ranging from 
190 yis to 1700 ps that trigger spark sources. Figure III-B-5 shows the 
delay unit circuitry. Figures III-B-6 and -7 shows typical shadowgraphs 
taken along the bullets’ trajectory before impact. The complete ballis- 
tic range configuartion is shown in Figures III-B-8 and -9* The bullet 
catcher is shown in Figure III- B- 10. 
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FIGURE m-B-l 

BALLISTIC RANGE COMPONENTS 



Figure III-B-' 



Ballistic Ran^e 



R-fle fyrte- 
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1 1 1 Ballistic Range Rifle Mount 







Figure III-B-4 Ballistic Range Shadowgraph Station. 
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V C c (5 Volts d.c.) 




FIGURE m-B-5 DELAYED PULSE GENERATOR 

CIRCUITRY (VARIABLE) 
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Figure III-B-6 Shadowgraph of a 63gr Semi-Point with a 
Velocity of 2941 fps . 




Figure III-B -7 Shadowgraph of a 50gr Spitzer Used in 
Range Calibration with a Velocity of 
3080 fps. 
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Figure lll-B-8 Ballistic Range Configuration (Down Range View) 





Figure III-B-9 



Ballistic 



Range Configuration 



I Up Range View). 
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Figure III-B-10 Ballistic Range aullet Catcher. 




C , Experimental Apparatus 



The experimental program was divided into three phases. The first 
phase consisted of measurements of the projectile energy loss during 
impact with the wall of an empty tank. Phase II measurements were taken 
with a water- filled tank but with the projectile entering the tank 
through a prepunched hole. Phase III measurements were taken using 
solid entry walls with a fluid- filled tank. A comparison of these 
results make it possible to partition the energy losses of each of the 
various hydraulic ram component dissipative processes. Figure III-C-1 
shows the test apparatus for Phase I. A test plate is positioned with 
two bullet sensors mounted directly behind. The average velocity after 
plate penetration is measured and the interaction phenomenum may be 
studied using a shadowgraph. 

Phase II and III testing uses the test setup shown in Figure III-C-2. 
Two shadowgraph stations provide for a photographic record of hydraulic 
ram phenomenum at two preset times after projectile impact. 

Figure III-C-3 shows the plexiglass side wall tank with its removable 
entry and exit walls. The entry and exit walls are fastened to the 
1 inch plexiglass side walls by studs that pass through steel clamping 
frames which sandwich the test plate. The tank is a cube with l8 inch 
inside dimensions. The electronic delay and spark source trigger cir- 
cuits used for the tank shadowgraph are similar to those shown in 
Figure III-B-5 . Station one timing starts at the last downrange shadow- 
graph station sensor and its delay circuit range is variable from 200 ^s 
to 1650 \is . Station two timing starts when the station one spark source 
fires and another delay circuit is used to fire the station two spark 
source from 3 to 20 jjls later. 
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Figure III-C-1 Test Specimen Setup. 







Station 1 




Schematic of Shadowgraph Apparatus 





Station 1 



Station 2 



Fig. III-C-2 Fuel Tank Test Setup 
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Fig. III-C-3 



Test Tank and Stand 
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IV. TEST RESULTS TO DATE 



A. Penetration Studies 

Phase I testing was conducted to gain some useful knowledge of the 
fuel tank wall-projectile penetration characteristics. Metal plates of 
7075- T6 aluminum of thicknesses o.05, 0.09, and 0.l6 inch were used. 
Projectiles of three masses and three nose shapes were fired into these 
plates at varying initial velocities using a .222 Remington rifle. The 
projectile velocity was measured before and after the plate penetration. 
All projectiles studied were stable before impact and entered the test 
plate perpendicular. Several shots at each mass, nose shape, and velocity 
were made to provide statistical data in order to determine the energy 
lost by the projectile during penetration of the entry wall without 
fluid support. The energy loss thus determined will be the same as the 
energy lost during a similar event in hydraulic ram except for the 
effects of fluid foundation and damping. Phase II testing will yield 
the amount of energy initially lost to the fluid with entry wall- fluid 
interactions but with no penetration effects. Phase III testing will 
yield the total initial energy loss during penetration. Comparisons of 
Phase II and Phase III should yield the energy loss during entry wall 
penetration with fluid effects. Comparison of these results with Phase 
I should then yield the difference in energy loss due to the presence 
of the fluid medium. 

To date the Phase I testing has been accomplished (Ref. 8) for the 
0.09 inch thickness plate only. At the present time test data for the 
other two wall thicknesses is being taken. Plans also include similar 
testing using 30 caliber projectiles in the near future. 

Figure IV-A-1 shows a typical projectile flowfield after plate 
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penetration. Spallation occurs from the back face of the entry wall 
producing small supersonic fragments that travel behind the projectile 
bow shock. A nearly spherical strong shock wave propagates into the 
tank from the projectile point of entry during the initial stage of plate 
penetration. Shortly after bullet passage through the plate, the bullet 
bow shock wave catches up to the original spherical wave. These waves 
coalesce and the resulting bullet bow shock wave has a reflex shape. 

Previous work by Forman (Ref. 9) derived equations for the minimum 
penetration velocity (ballistic limit) and the projectile exit velocity. 
The ballistic limit is given by: 



where G is the material shear modulus, r is the radius of the hole 
produced in the plate (assumed to be the same as the radius of the pro- 
jectile), t the wall thickness, and m is the projectile mass. For impact 



velocity. Assuming a constant energy loss the exit velocity may be 
found by: 



of nose shape and assumed that the kinetic energy loss is independent of 
impact velocity for impact velocities greater than the ballistic limit. 
Figure IV- A- 2 shows the results of firing projectiles of constant mass 
but different nose shapes. The data for both shapes agree closely with 
the theoretical predictions. This indicates that nose shape has a minor 
role in determining the exit velocity. Figures IV-A-3 and -4 show the 




( IV-A-l) 



velocities greater than V q the exit velocity V g approaches the impact 




(IV-A-2) 



where V_^ is the impact velocity. Equations (l) and (2) are independent 
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result of tests using projectiles of different mass with the same nose 
shape. These tests again showed excellent agreement with theory over 
the tested range of impact velocities. The greater length of the 63 gr. 
projectiles increase the effect any yaw or spin axis nutation has at 
plate impact. If entry is not perpendicular, damage is more severe and 
energy loss due to penetration is greater. A plot of kinetic energy loss 
versus initial velocity for all projectiles is shown in Figure IV-A-5. 

The theoretical energy loss is invariant with initial velocity and is 
shown as a reference. The energy loss associated with high mass project- 
iles is typically greater than that of the lower mass projectile for 
equal initial velocities as expected. 

Inspection of the damage to the entry walls shows that maximum 
damage occurs at a velocity slightly higher than the ballistic limit 
and the damage size approaches the diameter of the bullet as the impact 
velocity increases from the ballistic limit. 
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Figure IV-A-1 Flow Field after Plate Penetration. 





IMPACT VELOCITY, Vj (FPS) 



Fig, IV-A-2 Lx It Velocity versus Impact Velocity for Two 
22.2 Caliber Projectile Shapes Fired Normal 
to a 0.090 inch Thick 7075-T6 Aluminium plate. 
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EXIT VELOCITY, V p (FPS) 




I'ig. IV-A-3 Exit Velocity versus Impact Velocity for a 
55gr Semi-Pointed 22.2 Caliber Projectile 
fired Normal to a 0.090 inch Thick 7075-T6 
Aluminium Plate. 
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IMPACT VELOCITY, V; (FPS) 



Fig. IV-A-4 Fxit Velocity versus Impact Velocity for a 
63gr Semi-Pointed 22.2 Caliber Projectile 
Fired Normal to 0.090 inch Thick 7075-T6 
Aluminium Plate. 
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2000 2500 3000 

INITIAL VELOCITY, V; (FPS) 

Fig. iv-A-5 Plot of Kinetic Energy Loss versus Initial 
Velocity for a Range of Projectile Masses 
and Shapes. 



B. Shock Phase Studies 



Phase II and III studies of the shock phase of hydraulic ram to 
date have been limited to low energy .22 caliber shots. An attempt to 
increase kinetic energy level by using higher velocity .222 ammunition 
resulted in catastropic failure of the test tank as shown in Figure 
IV-B-1. The failure of this tank design to withstand hydraulic ram 
pressure loadings is a regrettably graphic example of the magnitude of 
the survivable fuel tank design problem. Inspection of the wall frac- 
ture shape seems to indicate that hydraulic ram tensile loads, concen- 
trated locally by the side wall stud threads, produced failure. It can 
also be deduced that since the front wall sustained the major damage, 
shock phase pressures can significantly reduce the ability of fuel tank 
front walls to survive the later phases of hydraulic ram. 

In spite of the test tank destruction, much useful information has 
been gathered on the shock phase of hydraulic ram. Dual spark shadow- 
graphs of the tank were taken at varying times after projectile impact. 
These data provided information on the characteristics of the shock 
wave pressure pulse and the rate of energy loss of the projectile. 

1. Shadowgraphs of Fuel Tank Flow Field 

Phase II testing was conducted using a 1 inch diameter pre- 
punched hole in a .09 inch thick entry wall. Figure IV- B- 2 illustrates 
typical results of the impact by a 1.95 gram lead bullet at a velocity 
of 0.374 kilometers per second (1226 ft/ sec) into the water- filled test 
tank. The photographs of Figure IV- B- 2 show that a hemispherical wave 
front is produced whose center is the point of impact. It is implied 
from this fact that projectile shape and forward momentum have no 
significant effect on the wave front motion during the early stages of 
hydraulic ram. Inspection of the wave thickness also shows that the 
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shock front strength varies along the wave. The greatest strength 
occurs in the direction of projectile flight and the least strength 
occurs where the shock and wall meet . 

Phase II testing involved impacts with standard factory loaded 
. 22LR ammunition. Average velocities of the resulting shock waves were 
determined from shadowgraphs similar to those in Figure IV-B-2. The 
shadowgraphs were used to determine time and position of the shock wave 
as it propagated through the fluid. The slopes of the resulting time- 
distance plots yielded propagation velocities of the shock. Figure IV- 
B-3 in such a plot of shock front velocity for a projectile with an 
averaged impact kinetic energy of 135.2 Joules. The plot shows that the 
shock front velocity decays to the acoustic speed in water (1.5 kilo- 
meters per second or 4920 ft per second) at a distance of about 11 
centimeters (4.3 in). 

Pressures generated in the water at the shock front corresponding 
to the measured shock front velocities were determined from data given 
in Table 2.2 in Ref. 10 . Figure IV-B-4 shows these shock front pressures. 
Extremely high pressures are generated in the water near the impact point 
and decay rapidly due to geometric expansion of the shock front. For 
example, Figure IV-B-4 shows that at a distance of 2.24 centimeters 
(0.88 in) from the impact point the resulting shock front pressure is 
3750 kilograms per square centimeters (53,335 lb/ ft ). As with any 
object moving through a fluid, the projectile creates a dynamic pressure 
in the water due to the stagnation conditions on the front surface. 

4 2 

Pressures of 710 kilograms per square centimeter ( 1.01 x 10 lb/ in ) 
were found to exist in the vicinity of the projectile during these 
studies. Figure IV-B-5 shows the evolution of the pressure field 
behind the shock front as a function of time. When the projectile 
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first penetrates the wall a strong shock is created and decays in 
strength as it propagates outward. At time t^ the pressure behind the 
shock greatly exceeds the stagnation pressure on the projectile nose. 
Comparing pressures at t^ and t^ it is apparent that both shock front 
pressure and projectile stagnation pressure decay; however, the decay 
of shock front pressure is much more rapid. At time t^ the stagnation 
pressure exceeds the shock front pressure. One could define the time 
where these pressures are equal as a characteristic time separating the 
shock and cavity phases of hydraulic ram. The characteristic time for 
the conditions of the test was found to be 32 microseconds. 

Figure IV-B-6 illustrates the later stages of cavity growth and 
projectile speed decay. The air- filled cavity behind the projectile 
is opaque to light and casts a sharp shadow. The true shape of the 
cavity is not that of the shadow because of light ray refraction through 
the test setup. The spark shadow of the cavity is wider than the true 
cavity size which could be found using X-ray shadowgraph techniques. 

The projectile shadow is not visible but its position may be inferred 
by the position of the bright cusp of light present at the leading edge 
of the cavity shadow. This cusp is formed by light bent around the 
projectile. This bending is caused by changes in fluid refractive index. 
It is expected that the position of the nose of the projectile corres- 
ponds to the tip of the cusp. Using this technique, measurements were 
made of the progress of the projectile through the tank after impact. 

A detailed analysis of projectile velocity decay is extremely diffi- 
cult because of the many variable parameters. Assumptions must be made 
to account for the changes in projectile shape, drag coefficient, and 
attitude. Because of these difficulties a simplified analysis was made. 
If the projectile drag is given by 
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(iV-B-l) 



the assumption of constant area, drag coefficient equal to one and 
constant projectile mass results in the following projectile velocity 
decay: 



Figure IV- B- 7 and IV- B- 8 show a comparison of experimental data with 
the previous analytical results for velocity and energy decay as a 
function of distance traveled by the projectile. This comparison indi- 
cates that in spite of the simplifying assumptions, the agreement 
between calculated and experimental values of velocity and energy after 
impact was good. In the experimental tests the projectile underwent 
a continual deformation during deceleration into the tank. The increase 
in frontal area during the deceleration was approximated by using an 
average final projectile frontal area of 0.5 square centimeters. 

Decay of the projectile kinetic energy can be directly related to 
the amount of energy being transferred to the fluid. The slope of the 
projectile energy curve gives the rate of energy transfer to the fluid, 
assuming no projectile deformation. The average energy transfer rate 
for the projectile tested was approximately 20.6 Joules per centimeter 
of fluid traversed. 

Phase III testing was accomplished using techniques similar to those 
used in Phase II. The test entry wall was however, solid. Figure IV- 
B-9 illustrates typical results of the impact by a 1.95 gram lead 




( IV-B-2 ) 
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bullet at a velocity of .37*+ kilometers per second (1226 ft/sec). The 
photographs of Figure IV-B-9 show the obscuration of events interior 
to the shock wave. The shadowgraphs show a strong hemispherical wave 
front as in the "pre-punched" wall case. This hemispherical wave is 
however, interacting with another wave front which appears as a straight 
line oblique to the wall. This wave front was caused by the propagation 
of shear and dilatation waves moving through the metal wall at a faster 
rate than the waves through the water. This wave is created by the 
penetration stage punching action of the projectile and its rate of 
propagation was found to be proportional to the bulk sonic velocity 
(dilatation wave) for aluminum (6.2 km/ sec.) 

2. Deformation of Lead Projectile 

Projectiles were recovered after each experimental test and 
were found to display a unique pattern of deformation. The deformation 
phenomenon was very consistent for projectiles impacting the fluid 
through pre-punched plates. Three examples of projectiles after this 
type of impact are shown in Figure IV-B-10. The extent of deformation 
was directly proportional to the impact velocity and all projectiles 
exhibited one flat side and an undeformed side. This type of deforma- 
tion is the result of the dynamic stability of the projectile as it 
traverses the tank fluid. As shown in Figure IV-B-11 the dynamic 
stability of the projectile is maintained through the equilibrium of 
the pressure and inertia. The figure demonstrates the effect of an 
angular displacement of the projectile from the equilibrium position. 

The center of pressure moves upward on a slightly curved front surface 
with the stagnation point. This creates a restoring moment which returns 
the projectile to its original position. Projectiles which penetrated 
a solid entry wall did not demonstrate this deformation phenomenon. 
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Figure IV-B-12 shows these projectiles and the small aluminum cap that 
was punched out by the projectile passage through the entry wall. 
Figure IV-B-13 shows typical entry wall damage. 
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Figure i V - H- 1 : 



Test Tank After Impact of .222 alit r 
High Energy Projectile 
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Figure IV-B-2: Shock Waves Resulting from a 

Prepunched Hole Impact 
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SHOCK VELOCITY (km/sec) 




Figure IV-B-3 Shock Front Velocity Generated in Water 
by Projectile Impact. 
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Figure IV-B-4 Shock Front Pressure Generated in 
Water by Projectile Impact. 
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Figure IV-B-5 Evolution of Pressure in the Transition from Shock to Cavity Phase. 
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Figure IV-B-6: Cavity Growth Resulting from a 

Prepunched Hole Impact 

66 




X ( cm ) 

FIGURE 3E-B-7 PROJECTILE VELOCITY DECAY RATIO 

AS A FUNCTION OF DISTANCE INTO 
TANK 
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FIGURE 33Z>B-8 EXPERIMENTAL AND CALCULATED 

RATIOS OF PROJECTILE ENERGY 
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Figure IV-B-9 Shock Waves Resulting from 0.22 Caliber 
Camber Projectile Impact 
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Figure IV-B-10 Projectiles After Impact Through 
Pre-punched Entry Wall 
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(B) Displaced from equilibrium 



Figure IV-B-ll Dynamic Stability of a Deformed 
Projectile in Water After Impact. 
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Figure IV-B-12 Projectiles After Impact Through 
Solid Entry Wall. 
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Figure IV-B-13 Solid Entrv Wall After Impact 
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V. RECOMMENDATIONS FOR FURTHER STUDY 



In FY 74, Hydraulic Ram Studies will be continued in the Department 
of Aeronautics. The following items are recommended as an outline for 
this future work: 

1. Shock phase shadowgraph studies should be extended to higher 
energy and caliber projectiles. Toward this end an experimental 
test tank capable of withstanding the damage caused by these 
projectiles has been designed. Information produced by the above 
studies should be correlated with analytical predictions of 
Yurkovitch. 

2. The test tank should be instrumented to measure internal pressures 
during the shock and cavity phases. This data will allow assess- 
ment of the Yurkovitch predictions during shock phase and - 
Lundstrom predictions during the cavity phase . 

3. The test tank entry wall should be instrumented with strain 
gages to measure wall stresses during hydraulic ram. This data 
will allow assessment of structural interaction prediction 
methods . 

4. The test tank entry wall should be instrumented to obtain wall 
motion during hydraulic ram. A fiber optic proximity gage is 
suggested as the measuring technique. This information will 
allow assessment of the structural interaction prediction method. 

5. Continuing studies should be made to ascertain the validity of the 
combined piston theory - fluid pressure analytical- numerical 
approach and to improve the methods where ever possible. 
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